Entanglement, avoided crossings and quantum chaos in an Ising 
model with a tilted magnetic field 
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We study a one-dimensional Ising model with a magnetic field and show that tilting 
the field induces a transition to quantum chaos. We explore the stationary states of 
this Hamiltonian to show the intimate connection between entanglement and avoided 
crossings. In general entanglement gets exchanged between the states undergoing 
an avoided crossing with an overall enhancement of multipartite entanglement at 
the closest point of approach, simultaneously accompanied by diminishing two-body 
entanglement as measured by concurrence. We find that both for stationary as well as 
nonstationary states, nonintegrability leads to a destruction of two-body correlations 
and distributes entanglement more globally. 

PACS numbers: 03.67.Mn,05.45.Mt 



INTRODUCTION 



ving states of various spin 



me evolving states 



10 



11 



12 



13| 



Entanglement content in the ground, thermal and^time^eyo 
models have been an active subject of recent research 
While many of the studies have been on integrable spin-chains such as the transverse Ising 
and the Heisenber g m odels, few have explored implications of nonintegrability and quantum 



chaos, for e.g. 
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ll| . Entanglement itself in these many-body systems is not characterized 



by a single number, but rather is revealed in various measures of multipartite correlations. 
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Motivated by the rapid developments in quantum information theory [14| these studies have 
revealed a rich phenomenology including entanglement scaling at zero temperature second 
order phase transitions and logarithmic (in number of spins) divergence of subsystem entropy 
at quantum critical points |s[ 6]. 

The relationship between quantum chaos and entanglement is a c omp lex one and has also 



been studied in various systems 
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271 ] . It appears that 



typically chaos can enhance entanglement especially of a multipartite kind. Results suggest 
that initially unentangled states are capable of developing large multipartite entanglement 
under quantum chaotic evolutions that are persistent in time 10, [20]. However integrable 
evolutions can generate large entanglements for specific initial states at specific times flol |. 
Early studies were based on bipartite systems such as coupled tops and pendula where 



numerical results 
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16j and random matrix modelling showed that the entanglement as 



measured by a subsystem von Neumann entropy was enhanced in regions of quantum chaos 



|16| . Later studies of many -body systems including spin models and quantum maps showed 
a subtler relationship 10, QjJ, |23|, |24j, but there is significant evidence that multipartite 



entanglement is enhanced by quantum chaos jlO|, [20|, [24|. Related systems are disordered 
spin chains wherein transitions in level statistics with increasing disorder have been shown 
to be correlated with decreasing two-body entanglements, for instance in 25]. 

The mechanisms that underlie the correlations between entanglement and nonintegrabil- 
ity are not entirely explored. It may simply be that the entanglement content of random 
states are reflected in systems with quantum chaos. The critical requirement of operators 
for producing large multipartite entanglement has been explored before 24| and there are 
indications that the random distribution of matrix elements, rather than other character- 
istics such as the nearest-neighbor spacing distribution (NNSD), are needed for generating 
high entanglement. On the other hand such a distribution may not be easily generated. 
Also connections between localization and entanglement have been noticed earlier 26l. 127]. 

In this paper we point to connections between nonintegrability and entanglement via 
avoided crossings. One of the hallmarks of quantum chaos is level repulsion, the tendency of 
quantum energy levels to avoid each other {28J. This is most clearly seen when a parameter 
of the system is changed. The resultant energy level dynamics has a typical behavior for 
quantum chaotic systems once the levels are restricted to the same symmetry class. There 



are no exact or accidental crossings, levels come close to each other and get "scattered". 
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The effect of this is reflected in the most popular diagnostic of quantum chaos, namely the 
nearest-neighbor spacing statistic. Due to level repulsion this deviates from the Poisson 
statistic^ obtained for integrable systems to the Wigner distribution typical of quantum 
chaos 



When two levels approach each other due to the variation of one parameter, from a theo- 
rem of von Neumann and Wigner we know that generically they will not become degenerate 



29j . The avoided crossing can be very sharp or broad and may also involve other levels. The 



nature of the eigenstates "morph" into one another at an avoided crossing 30]. As avoided 
crossings are generically found in nonintegrable systems, it seems natural to look at the 
behavior of entanglement at these points. This has been exploited to some extent in earlier 
works that seek to create entangled states such as the W or GHZ state by using the super- 



331 ] . we study this more systematically 



positions that develop at avoided crossings [31, 32], 
in a nonintegrable system. For this same system we show that there is a close correlation 
of the NNSD to the extent of entanglement; the Poisson statistic favoring nearest neighbor 
correlations and low multipartite ones, while the Wigner distribution favors large multipar- 
tite entanglements with low two-body correlations. We show that for time evolving states 
this is reflected in the way an initially maximally entangled pair of spins evolve and share 
this entanglement along with that created by the dynamics. Once again the case when there 
is chaos leads to a destruction of two-body correlations alongwith enhanced multipartite 
entanglement. 

II. THE MODEL 

The Hamiltonian we will use in this study is 

H(J, B, 9) = J J2 «+i + B £(sin(0K + cos(0)<) (1) 

n=l n=l 

For 9 = the magnetic field is longitudinal, the model is almost trivially integrable and 
the spectrum is highly degenerate. When 9 = 90° the field is transverse and the model 



is still integrable thanks to the Jordan- Wigner transform [34], |35( that maps the model to 
one of noninteracting fermions. This has been extensively studied both from the original 
statistical physics motivation and also recently from the viewpoints of quantum information 
theory. When < 9 < 90°, the model is not integrable, when converted using the Jordan- 
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Wigner transform the resulting fermions are interacting. When the magnetic field is pulsed 
or kicked, this has been studied as the "kicked Ising model" and has also been used to 
study entang.e.ent Q. Tbe MCd Ising m „de. is believed to be q _ chaotic for 
intermediate tilt angles j^]. However the time dependence is not essential in introducing 
nonintegrability and we will study the autonomous Hamiltonian above. It is evidently a 



very natural generalization of the well studied case of the transverse Ising model 35|. We 



have since doing this work noticed that it has appeared in two other recent complimentary 
studies [ll|, 3^- 

The antiferromagnetic and ferromagnetic chain spectrums are related to each other as 

(n®°?) H ^ B ^) (n®tf) = -H(-j,B,e). (2) 

The ground state of the antiferromagnetic chain is the spin-flipped version of the most 
excited state of the corresponding ferromagnetic chain. We will consider J > 0, as we will 
mostly deal with the entire spectrum of states. A discrete symmetry present for all angles 
is that of interchanging the spins at the sites % and L — i + 1 for all i = 1, . . . L, a "bit- 
reversal" symmetry as the field and interaction do not distinguish the spins except for the 
open boundary condition. If B represents this reversal 

B\s 1 82...8 N ) = \8 N ...8 2 S 1 ), [H,B]=0 (3) 

where \si) is any single particle basis state, such as the standard (S z ) basis states we will 
use. Since B 2 = 1, the eigenstates can be classified as odd or even on bit-reversal. The 
dimensionality of the odd subspace is half the number of non-palindromic binary words of 
length L, while the even subspace is larger by the number of palindromic words of the same 
length. The chain with periodic boundary conditions naturally has a much larger symmetry 
group corresponding to a shift operation. For 6 = the spectrum is highly degenerate 
with eigenstates that can be chosen to be product states. For the transverse Ising model, 
9 = 7r/2, the spectrum is much less degenerate, but the Jordan- Wigner transform converts 
the problem into one of non-interacting fermions. For intermediate angles this transform 
leads to a model of interacting fermions which is not solvable. 

In Fig. (CTJ we show a part of the even subspace of the energy spectrum for L = 8 spins. 
The fanning out of the energy eigenvalues from degenerate ones is typical of systems that 
lose symmetry and becomes nonintegrable. The apparent crossings of energy levels are in 
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FIG. 1: A part of the even subspace of the energy spectrum for the case of L = 8 spins. J = B = 1. 

fact very close avoided crossings. As the tilt angle increases further the avoided crossings 
become more apparent and dense. 

A standard unfolding of the spectrum in a given symmetry class is done by numerically 
fitting the staircase function (cumulative density) to high-order polynomials and using this to 
map the energies to unfolded ones such that the mean energy spacing is unity. The NNSD is 
a commonly used indicator of quantum chaos[28, 43]. In Fig. ([2]) we show the NNSD for four 
angles, three of them fairly close to a purely transversal field {0/tt = 99/200, 15/32, 7/16). 
It is to be noted that even for small longitudinal fields the NNSD is close to the Wigner 
distribution, corresponding to the Gaussian Orthogonal Ensemble (GOE) of Random Matrix 
Theory (RMT) |28|,|43j, given by 

7T 



Pw(s) 



s exp( — its /4) 



The solvable pure transverse field case is not apparently completely desymmetrized by 
the "bit-reversal" operator, degeneracies persist and therefore we do not plot the NNSD in 
this case. As the tilt angle is increased further we are tending towards the solvable purely 
longitudinal case, with very large degeneracies. A reflection of this is seen in the NNSD 
deviating significantly from the GOE distribution even at 9 = tt/3. The appearance of the 
GOE distribution indicates a certain measure of quantum chaos and applicability of random 
matrix ensembles. However it must be noted that the Hamiltonian we have considered is 
"simple" in that there is no disorder and the interaction is only nearest neighbor. We will 
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FIG. 2: The nearest neighbor spacing distribution of the even states of a chain with L = 13 spins 
at various angles 6 of tilt of the magnetic field. J = B = 1. 

presently quantify the distance of the NNSD from the Poisson one (Pp = exp (— s)) and that 
of the GOE (P w (s)). 



III. ENTANGLEMENT IN THE STATIONARY STATES 

We will use three measures to quantify entanglement. 

1. Entanglement within pure states of a bipartite system can be measured by the von 
Neumann entropy of the reduced density matrices. However there are many bipartite 
splits, given L spins. Generalized entanglement measures have been constructed and 
studied based on the entanglement in each of these possible splits. We will however 
consider only the entanglement of the first L/2 spins with the remaining L/2 ones for 
L even, denoted as Sl/2- Thus if \ip) is an L— spin pure state, 

S L /2 = -tri,..., L/2 (pi,...,l/2 log(pi,...,L/ 2 )) (4) 

where 

Pi,..., L/2 = tr L/2+ i,...,L (I^X^I) (5) 

2. The state of any two spins, such as nearest neighbors, is in general a mixed state. For 
such a state while the entanglement can be measured as the average entanglement of its 



7 



pure-state decompositions, the existence of an infinite number of such decompositions 
makes their minimization over this set a nontrivial task. Wootters and Hill [38| carried 
out such a procedure for the case of two spin one-half (qubit) systems and showed that 
a new quantity they called concurrence was a measure of entanglement. This facilitated 
the study of entanglement sharing among many qubits. If pij is the reduced density 
matrix obtained by tracing out all spins except those at sites % and j, and defining a 
spin-flip operator, which takes p^ to 

Pij = K ® a v)Pij K ® a v) > ( 6 ) 
the concurrence of p i:; - is defined to be: 

C(pij) = max - \f\~2 - \fh ~ \[>4, 0} (7) 

where A, are the eigenvalues of the non-Hermitian matrix PijP~ij. Wootters {ssl showed 
that the entanglement of formation of pij is a monotonic function of its concurrence 
and that as the concurrence varies over its possible range [0, 1], the entanglement of 
formation also varies from to 1, We will use the square of the concurrence, called 
tangle, summed over all possible pairs of spins 

(r) = E C\ Pt3 ) (8) 

i<j 

as a measure of two-body correlations in the L— spin state. 

3. The Meyer and Wallach Q measure is also being widely used as a measure of multipar- 
tite entanglement. The geometric multipartite entanglement measure Q (3^], has been 



shown to be simply related to one-qubit purities [40], which makes their calculation 
and interpretation straightforward. If pk is the reduced densty matrix of the spin k 
obtained by tracing out the rest of the spins then 

L 



(9) 



In Fig. ([3]) we show the entanglement content in eigenstates across a particular spec- 
trum. We see that there is a secular dependence of the multipartite entanglement on the 
energy. We connect this dependence with other more grosser features such as localiza- 
tion of the eigenstates in the computational basis. Thus if ipk = (k\ip) is the k-th eigen- 
state component, k = 0,...,2 L — 1, we measure its spread in the computational basis 
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FIG. 3: The log of the participation ratio (solid line, top) and 4Q (dashed line, bottom) versus 
energy are shown for all eigenstates of an integrable (0 = tt/2) and a nonintegrable (9 = 7it/1Q case 
in the top panels. In the bottom panels are shown, for the same cases, the entanglement measure 
Sl/2 versus the Shannon entropy S s h of the same states. Notice the apparent lack of correlation for 
the integrable case as compared with the nonintegrable one. In all cases L = 10 and J = B = 1. 

as log(P-R) = -log(£* |^fc| 4 ) or as S sh = -Efcl^fc| 2 log(|^fc| 2 ), i.e., the log of the par- 
ticipation ratio and the Shannon entropy of the states. Both these measures are obviously 
single-particle basis dependent while measures of entanglement are immune to these changes. 
However there is a strong correlation between measures of entanglement such as Q and these 
grosser measures of localization, especially in the nonintegrable regimes. It is also clear from 
these graphs that average entanglement over the spectrum will be dominated by states from 
the central regions, away from the ground states of the antiferromagnetic chain under con- 
sideration, as well as away from the groundstate of the ferromagnetic chain to which the 
highest excited states correspond. 

We now turn to how individual states' entanglement content change as the parameter 9 is 
varied. Parametric motion of energy levels has been the subject of several studies in the past 
and their connection to quantum chaos is well known. To be specific we look at three energy 
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levels undergoing multiple avoided crossings in a range of the parameter, some of them 
sharp and some soft. Shown in Fig. (j3J) are the levels alongwith three measures of their 
entanglement, a sum of all two-tangles, the Q measure and the entropy Sl/2- The central 
state undergoes collisions with its neighbors, while the other two also undergo collisions 
with states not shown in the figure. It is clear from all the measures of entanglement that at 
avoided crossings there is an exchange of entanglement. In particular it is interesting that 
the highest energy state shown has a vanishing entropy Sl/2 till it undergoes an AC with the 
middle level which subsequently has zero entanglement, which further gets transferred to 
the state with the lowest energy. Thus along with avoided crossings at which states undergo 
structural changes, it is reasonable that entanglement properties get exchanged as well. In 
the vicinity of an AC the average entanglement of the two states involved increases sharply 
in ways that seem to depend on how sharp the crossing is. However at an AC while the 
multipartite entanglement increases, the two-body correlations decrease. This is illustrated 
in Fig. ([5]), where the total 2-tangle is seen to dip to almost zero at the point of closest 
approach of the AC. We note that a fine sweeping of the parameter in the vicinity of the AC 
is necessary to pick up this trend which is consequently not obvious in a coarser one such 
as in Fig. (TJJ 

The fact that at an AC entanglement is enhanced, and that at the point of closest ap- 
proach the entanglement is maximized has already been used before in the context of few 
spins and the GHZ state. However we see it here in the context of level dynamics of a 
nonintegrable system how generic it is, and that there are tradeoffs between multipartite 
and two-body entanglements. It maybe one mechanism by which nonintegrability enhances 
multipartite entanglement. A more detailed study of entanglement at ACs is called for, 
especially considering that what happens at such points is a coherent superposition of or- 
thogonal states. That superpositions will always cause an increase of entanglement is of 
course untrue, however it could be true for some large class or measure of states that are 
relevant to ACs. A very recent study of entanglement due to superpositions has appeared 
4=11 ] . and maybe relevant to what happens in the vicinity of an AC. 

We add two other numerical results that show the effect of nonintegrability on entangle- 
ment. We calculate the entropy Si of the first / spins of the chain. This is known to scale as 
log(Z) for large I at quantum critical points and tends to a constant otherwise. However it is 
also known that random states of relevance to quantum chaotic systems have subsystem en- 
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FIG. 4: clockwise from top-left: (1) details of three energy levels from Fig. ([I]) (J = B = 1,L = 
8) undergoing avoided crossings, (2) The sum of the 2-tangles of all spin pairs, (3) The global 
entanglement measure Q for the three states and (4) the entanglement of one-half of the chain 
with the other {S^^). 

tropies that scale as log(D) where D is the dimensionality of the subspace. If we are dealing 
with collections of qubits then D = 2 l and such subsystem entropies should scale linearly 
with /, Indeed in the Ising model under consideration along with a transition to classical 
chaos we observe a transition to a linear behavior. We show this in Fig. ([6]) where we take 
data from one hundred states at the center of the spectrum to smooth out state-to-state 
fluctuations. Since we cannot numerically access very large spin chains we show calculations 
based on a chain of 14 spins. The entropy S L / 2 shown previously is a special case of Si, 
but also appears to deviate from the straight lines shown in these graphs, deviations most 
likely originating in the symmetry of the spin chain. As a more direct connection to quan- 
tum chaos we quantify the distance of the spectral NNS distributions from the Poisson and 
the GOE distribution of RMT using the Kolmogorov-Smirnov (KS) goodness-of-fit statistic 
and correlate this to entanglement. The KS statistic is defined as the maximum difference 
between a cumulative distribution of the data and that of a hypothesis. As Fig. (J7|) shows 
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FIG. 5: The average Q and average total 2-tangle across two avoided crossings shown in Fig. (j4]) 
between the lower two levels, one soft and one hard. 
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FIG. 6: The subsystem entropy of Si of the first I spins in a chain with L = 14 spins (J = B = 1) 
for various tilt angles. 
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FIG. 7: The subsystem entropy of S1/2 calculated as an average over all the eigenstates, in a chain 
with L = 12 spins ( J = B = 1) for a range of tilt angles of the field, the maximum corresponding 
to the transverse case. The inset shows the Kolmogorov-Smirnov goodness-of-fit statistic for the 
corresponding spectral NNS distributions, the top line corresponding to the test hypothesis being 
the Poisson distribution, while the lower one corresponds to the GOE distribution. 



as the field gets tilted away from the transverse case, there is an increase in the entangle- 
ment as measured by Sl/2 while simultaneously the NNSD makes a transition to the GOE. 
This is seen for other global measures of entanglement such as the Q measure as well. This 
quantifies the transition from Poisson to GOE distributions of the NNS as well as correlates 
with chain entanglement properties. A similar correlation is seen when the angle approaches 
zero and the there is once more a deviation from the GOE. 



IV. TIME EVOLVING STATES 



In this section we investigate the dynamics of entanglement for the above model. Our 
strategy is to start with an initial state with exactly one maximally entangled pair and the 
others in an unentangled state and study its entanglement properties under time evolution. 
The initial state \ip (0)) with given entanglement properties evolves to \ip (t)) with a com- 
plicated distribution of entanglement. We investigate how the entanglement distribution 
in the chain varies with time by studying the evolution of nearest neighbur entanglement 
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(calculated using the concurrence measure), and the multipartite entanglement using the Q 
measure in a chain of 10 qubit. We consider a chain of 10 qubits with B = 1 and J = 1 
(units chosen so that h — 1) and with the initial state: 



The above state evolves with time and the plots below depict the distribution of entanglement 
across the chain function of time. 

Fig. [8] shows the density plot of the nearest neighbor entanglement (calculated using 
concurrence) as a function of /, the location of the pair, and time, for a chain of 10 qubits 
with \ip e ) as the initial state. The darkness of a region is proportional to entanglement for 
the given pair. We observe that entanglement is initially concentrated only between the 
first two qubits but as time progresses it gets distributed among other nearest neighbor 
pairs in the chain. We see that there is a threshold time which occurs fairly early into 
the evolution at which all the nearest neighbor pairs acquire entanglement. In some other 
studies that have investigated such transport, such as for the Heisenberg XY spin chain, the 
entanglement transport has been a gradual phenomenon |42j and tends to move as a patch 
across the chain with a constant velocity before being dispersed across the chain. 

However, in this case there is a sharp increase in the entanglement of the nearest neigh- 
bors and persists for a while before disappearing. The presence of this feature across the 
values of 9 and number of spins leads us to believe that it is a robust feature of the Hamilto- 
nian considered. The dynamics of the Q measure plotted in Fig. [9] also reveals interesting 
features. The sudden generation of nearest neighbor entanglement initially is accompanied 
by a steep rise in the Q measure. Since the Q measure quantifies the multipartite entan- 
glement in the chain, it is very likely that this steep rise is due to the nearest neighbor 
entanglement alone. This is followed by a drop in the nearest neighbor entanglement while 
the Q measure remains at a large value, indicating that as time progresses the entanglement 
is exclusively multipartite in nature. Further proof is obtained from the plot of the average 
two-body entanglement versus time. This quantity falls with time showing that the pairwise 
correlations are reduced and therefore we conclude that almost all the entanglement that 
survives after the initial phase is purely multipartite in nature. 

We notice a close connection between the non-integrability introduced by 9 and the 
entanglement dynamics. Specifically, we notice that for those values of 9 that result in high 
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FIG. 8: The nearest neighbor concurrence is plotted with time for four different tilt angles {6 /it = 
1/2, 5/12, 1/3, 1/6, from left to right and top to bottom) in a chain of 10 qubits. On the y-axis the 
scale is 1 = 1/25 time units and J = B = 1. 

non-integrability for the system, the time taken for the quenching of two-body entanglement 
(concurrence) is considerably lesser. Here, the degree of non-integrability is measured by the 
closeness of its energy spacing distribution curve to the Wigner distribution. For 9 = 7r/3, 
the spacing distribution is very close to Wigner and we observe a quick death of two body 
correlations. Concurrence is seen to persist for a longer time in cases where the spacing 
distribution is significantly different from the Wigner distribution, for e.g. 9 = n/2 or 
7r/6. Further, we also noticed in calculations not presented here, that average tripartite 
entanglement characterized by negativity also falls down faster for systems with greater 
non-integrability though it persists for a somewhat longer time than the corresponding two 
body correlations. Therefore it seems likely that the presence of non-integrability favors the 
generation of higher order correlations at the expense of lower order ones. 

To summarize, we have studied an Ising spin chain with a tilted magnetic field of which 
the well studied one-dimensional transverse Ising model is a special case. We have shown an 
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FIG. 9: The Q Measure is plotted with time for the same four tilt angles 8 as in the previous figure 
for a chain of 10 qubits and J = B = 1. 

intimate relationship between the spectral fluctuations of the Hamiltonian and the entangle- 
ment properties of the eigenstates. In particular we have shown that there is a larger mul- 
tipartite entanglement content when the spectral fluctuations are described by the Wigner 
distribution usually found for nonintegrable systems. At the same time, such chains are low 
on the nearest neighbor two-body correlations such as concurrence. We have also shown that 
at avoided crossings when there is generically a superposition of "bare" states, multipartite 
entanglement is enhanced while the two body concurrences are lowered. The properties of 
stationary states are reflected in those of nonstationary states, and in particular we have seen 
that when the chain is nonintegrable, initially present two-body entanglements are quickly 
quenched and are converted to multipartite entanglement while for integrable chains, signif- 
icant entanglement is retained in the form of two-body concurrences. 
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